Percolation line of stable clusters in supercritical fluids 
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We predict that self-bound clusters of particles exist in the supercritical phase of simple fluids. 
These clusters, whose internal temperature is lower than the global temperature of the system, 
define a percolation line that starts at the critical point. This line should be physically observable. 
Possible experiments showing the validity of these predictions are discussed. 
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Recent studies |jj of the supercritical phase of fluids 
reveal new and unexpected phenomena that concern a 
variety of domains, ranging from fundamental problems 
in cluster physics M to industrial applications S. A 
question relevant to many of these domains is the possi- 
ble existence of clusters of particles. We predict in this 
Letter that self-bound clusters of particles exist in the 
supercritical phase of simple fluids and that these clus- 
ters define a physically observable percolation line that 
starts at the critical point. 

Clusters of particles in the supercritical phase have 
been considered in the past, mainly from a theoretical 
point of view However, these clusters are math- 

ematical objects, introduced to match the percolation 
threshold with the thermodynamical critical point. Fur- 
thermore, any critical percolation line defined by these 
clusters has been considered as unphysical and efforts 
have been made to eliminate them by a proper re- 
definition of the clusters [[U?]] . 

The results presented here are based on extensive 
numerical calculations || of the cluster size distribu- 
tions obtained from canonical Monte Carlo and micro- 
canonical Molecular Dynamics simulations of systems of 
TV < 11664 particles which interact through the Lennard- 
Jones potential: V(r,j) = 4e[(r /ry) 12 - (ro/r#) 6 ]. The 
cutoff distance of the potential is fixed at 3ro and parti- 
cles are confined in a cubic box with periodic boundary 
conditions. The phase diagram of this fluid is represented 
in figure 1 PJIofl. This phase diagram as well as the re- 
sults we will show below are insensitive to this particular 
choice of the potential and can be considered generic for 
simple fluids. 

We assume that the system is at equilibrium when 
mean potential (kinetic) energy reaches a stable value 
in the canonical (microcanonical) ensemble. In the latter 
we take 2/3 of the mean kinetic energy as a measure of 
the temperature of the system. 
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FIG. 1. Phase diagram of the Lennard- Jones fluid. The 
dashed line indicates the boundaries of the gas (G), liquid (L) 
and solid phases (S) /P|7fc|/. The continuous line (main pur- 
pose of this Letter) is the critical percolation line calculated 
with the condition of stability of clusters. Temperature T and 
density p are in units of e and r^~ 3 respectively. 

The cluster size distributions are calculated once ther- 
mal equilibrium is reached. Physical clusters are defined 
according to a prescription proposed by Hill [jll|: At a 
given time, two particles are linked if their potential en- 
ergy exceeds their relative kinetic energy. A set of linked 
particles forms a cluster. We have checked that using 
this procedure, for equilibrium configurations, most clus- 
ters are stable by particle emission, namely the kinetic 
energy of each particle relative to the center of mass of 
the cluster is less than the sum of the potential energies 
due to the other particles of the cluster. 
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FIG. 2. Cluster size distributions n(s) at points (a)-(b), 

(c) -(d) and (e) of the phase diagram of figure 1. For curve 

(d) , the contribution of the largest cluster S m ax is sharply 
peaked beyond the boundaries of the figure. 

The mean cluster size distributions n(s) (s = 
1,2- • -,7V) calculated at the points (a,b,c,d,e) of the 
phase diagram are displayed in figure 2. One remarks 
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that at points (a) and (b) the distributions can be fit- 
ted (for s > 3) by the same power law distribution 
n(s) ~ s~ T with t — 2.20 ± 0.05. In fact, the same 
power law behavior is found along the full line of figure 
1. This line T p (p) is defined as the locus of points where 
m 2 reaches its maximum value |l2] ] . The are the usual 
moments of n'(s) [Q, where n'(s) is the distribution of fi- 
nite size clusters, excluding the largest one. We find that 
to the right of this line there exists a percolating cluster. 
Its size S max /N vanishes as (T p (p) —T) 13 . In the vicinity 
of that line the moments behave as m 2 ~ |T — Tj,(p)| -7 
and m 3 /m 2 ~ \T - T p (p)\- 1 / a . The values we find for 
these exponents are in agreement with those of random 
percolation theory (see table 1.) A more rigorous de- 
termination of the ratios of critical exponents can be 
made using finite size scaling relations. As a function 
of the size N of the system, one expects just on the crit- 
ical line the behaviors: S max ~ N 1 ~' 3 / 3u , m 2 ~ N~</ 3u , 
m 3 /m 2 ~ N l / 3a ", < Sl ax > - < S max > 2 ~ N 1+ ^ 3 » 
fL3f , where v is the critical exponent associated to the 
two-point correlation function Q. These quantities, cor- 
responding to point (b) of the phase diagram, are dis- 
played as a function of N in figure 3. One observes that 
beyond N ~ 50 the expected power law behavior is very 
well satisfied. The values we find for the ratios of critical 
exponents are also in good agreement with those of ran- 
dom percolation (table 1.). We get the same agreement 
at other points of the line. These results strongly sug- 
gest that clusters of self-bound sets of particles define a 
critical percolation line ( CPL ), characterized by the uni- 
versal exponents of random percolation. In addition, the 
present numerical calculations show that, within the un- 
certainties inherent to the finite size of the system and 
to the critical slowing down, this percolation curve ends 
at the (thermodynamical) critical point (T c ,p c ). Just at 
this point, for the same technical reasons, we are un- 
able to calculate accurately the corresponding critical 
exponents. Extrapolating the exact results of Kastclcyn- 
Fortuin (lj] and Coniglio-Klein for the lattice gas model 
0, one expects to get at the (thermodynamical) critical 
point (and only at this point) the critical exponents of 
the Ising model. 
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Present Work 
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0.41 


0.4 ± .1 
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1.80 


1.6 ±.3 


a 


0.45 


0.44 ± .05 


T 


2.18 


2.20 ± .05 


0/v 


0.47 


0.48 ± .02 


j/v 


2.05 


2.0 ±.05 


1/au 


2.53 


2.60 ± .05 



TABLE I. Critical exponents associated with the cluster 
size distributions, in 3d Random Percolation jfy and present 
work. 



At p < p c or T < T c we do not see any signal of critical 
(percolation) behavior. The distribution n(s) decreases 
much faster than s -2 - 2 and clearly deviates from a power 
law form. However, a "macroscopic" cluster appears as 
soon as one penetrates into the two-phase region. This is 
clearly seen in figure 2, where we plot two distributions 
n(s) around point (e) of the phase diagram, with temper- 
atures of T = 1.10 and T — 1.05, just above and below 
the liquid-gas coexistence curve. This sharp signal corre- 
sponds very well to the crossing of the coexistence curve. 
Indeed, small change in T induces a drastic change in 
n(s) for large s |l5| . 

We have checked that all the above findings do not de- 
pend on the specific definition of stable clusters: Using 
two different definitions, Jl6|,[l7j , both based on minimiza- 
tion procedures of the interaction energy between differ- 
ent clusters, we get cluster distributions that are almost 
identical @. 

It is interesting to emphasize that the CPL, character- 
ized by the critical exponents of random percolation, is 
found, on a deterministic dynamical framework, without 
any explicit reference to a random (site or bond) perco- 
lation mechanism. This is even more striking when using 
definitions [ p^JT^ ] , for which clusters result from a global 
energy balance (and not simply from a bond activation 
prescription). 




10 100 1000 10000 

N 

FIG. 3. Finite size scaling behavior of Smax, mi, m^/mi 
and a 2 =< Smax 2 > — < Smax > 2 calculated at point (b) of 
the phase diagram. The straight lines are the best fits in the 
range 50 < TV < 11664. 

A this point, a question that arises naturally is the 
link between this CPL and the critical percolation line of 
references . We recall that the latter is found in the 
lattice gas model when choosing the Coniglio and Klein 
definition of clusters j| . In ref. |l^] it was shown that this 
definition is (nearly) equivalent, in a lattice gas model, to 
a condition of stability of clusters by monomer emission. 
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In this sense, the CPL can be seen as a generalization to 
a realistic fluid. 

It is also interesting to recall that, when defining clus- 
ters with a real space criteria, for instance as sets of par- 
ticles that are two by two at a distance less than some 
distance d c , a percolation line also exists p|, pOpl| ], but 
its position depends crucially on the choice of d c . How- 
ever, this percolation line may be relevant in situations 
in which physical clusters are sets of particles close in real 
space, for example in conductivity experiments in which 
neighboring atoms exchange electrons in overlapping or- 
bits. 

As a natural consequence of the definition of clusters 
based on a stability condition, we find that the internal 
" effective temperature" T e ff of clusters Q is always 
less than the global temperature T of the system. Using 
the molecular dynamics calculation, we find that T e ff 
grows as a function of s from zero for s = 1, to a limiting 
value for s > 100. In the supercritical phase, this limiting 
temperature, that depends essentially on the density of 
the system, ranges between T e // = 0.7 and T e ff = 1.0, 
i.e. small clusters behave as "solid-like crystals" while the 
large ones behave as "liquid-like droplets" , as can be seen 
considering the complete phase diagram. We have also 
observed that along the CPL the total energy of the sys- 
tem remains almost constant. The origin of this energy 
invariance, that results from a subtle balance between 
internal, center of mass kinetic energies and inter-clustcr 
potential energies, is not understood. 

Finally, we briefly discuss the possibilities to observe 
experimentally the CPL. Notice first that this line should 
not be confused with the Fisher- Widom line [^3| , or with 
the extrapolation of the rectilinear diameters line pi[ . 
Those two lines divide the supercritical region into a gas 
like and a liquid like domain with probably no relevance 
on clustering. 

Many experiments dealing with critical behavior of bi- 
nary fluids (which belong to the same universality class as 
liquid-gas and 3d Ising model) have been performed [^5) . 
Observations of concentration fluctuations in the mix- 
ture of isobutyric acid and water have been done in the 
vicinity of the critical point |26]| and a fractal dimension 
Df = 2.8±0.1 of clusters has been determined. The value 
we determine along the CPL is Df = \jav = 2.60 ±0.05. 
However, the relationship between these clusters, defined 
from persistent density fluctuations and the stable clus- 
ters is not yet clear. Work is in progress in this direction. 
The mobility of H + ions in the same binary mixture has 
also been studied very recently Effl. The sharp decrease 
of the ion mobility observed for T > T c as the critical 
concentration is exceeded, has been associated to the ap- 
pearance of a percolation line of dynamical clusters. The 
position of this curve in the p — T plane would suggest 
however a closer connection with clustering in real space 

In a different domain, early experiments on the effusion 



of a fluid through a pin hole [|28| have shown the presence 
of stable clusters (dimers and trimers) in the gas phase. 
In these experiments, the mass yield of escaped clusters 
with a given velocity can be related to the moments rrik 
of the cluster size distribution n'(s). Unfortunately, the 
interesting regime, namely the crossing of the CPL has 
not yet been explored. It should also be possible to iso- 
late stable clusters by a sudden disassembly of a piece 
of fluid into droplets, like in the fragmentation of atomic 
nuclei or atomic aggregates |2^,|3(| . Several experiments 
p2| , p9| show indeed, as a function of the excitation en- 
ergy, an evolution of the fragment size distributions n(s) 
that suggests a crossing of the CPL. However, due to the 
small size of the system, the results are still inconclusive. 

In summary: Large scale and long time Monte Carlo 
and Molecular Dynamics calculations suggest that en- 
ergetically stable clusters are present in a supercritical 
simple fluid. These clusters, that are cooler than the 
system, percolate below a critical line that ends at the 
(thcrmodynamical) critical point. Various experiments, 
within the reach of present day experimental techniques, 
are suggested to show the presence of these stable clusters 
and the existence of this critical percolation line. 

We have benefited from fruitful discussions with O. 
Martin and E. Plagnol. O. Bohigas is specially acknowl- 
edged for his critical remarks. 
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